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The carrier spin and impurity spin densities in diluted magnetic semiconductors are considered 
using a semiclassical approach. Equations of motions for the spin densities and the carrier spin 
current density in the paramagnetic phase are derived, exhibiting their coupled diffusive dynamics. 
The dynamical spin susceptibilities are obtained from these equations. The theory holds for p-type 
and n-type semiconductors doped with magnetic ions of arbitrary spin quantum number. Spin-orbit 
coupling in the valence band is shown to lead to anisotropic spin diffusion and to a suppression of 
the Curie temperature in p-type materials. As an application we derive the Hall-voltage noise in 
the paramagnetic phase. This quantity is critically enhanced close to the Curie temperature due to 
the contribution from the anomalous Ifall effect. 

PACS numbers: 75.50.Pp, 75.40.Gb, 72.20.My, 03.65.Sq 
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I. INTRODUCTION 

In recent years, a lot of progress has been made in 
the physics of diluted magnetic semiconductors (DMS), 
in particular in III-V materials doped with manganese. 
In the best studied material, (Ga,Mn)As, ferromag- 
netic transition temperatures around 160 K have been 
achieved. ^-^ On the theoretical side, a Zener model based 
on valence-band holes exchange-coupled to local impu- 
rity spins is very successful in describing this material, at 
least in the metallic regime i^i^i^iSii In (Ga,Mn)As man- 
ganese acts as an acceptor and introduces localized spins 
S = 5/2 due to its half- filled d-shell. The material is p- 
type but partly compensated, probably due to arsenic 
antisitesS^ and manganese interstitials.''^''' In group-IV 
semiconductorsii manganese plays a similar role. On 
the other hand, in II- VI materials manganese introduces 
a spin but is isovalent with the host cations. 

It has also been realized that disorder is crucial for 
the understanding of the properties of DMS, even in the 
metallic regime liSiiiiiii^ There are two main scattering 
mechanisms: disorder scattering due to the Coulomb 
potential of charged donors and acceptors and spin- 
exchange scattering off randomly distributed impurity 
spins. The Coulomb interaction is the dominant contri- 
bution to disorder. This is due to compensation, which 
leads to a lower hole concentration and thus on the one 
hand to the presence of charged defects of either sign and 
on the other to less effective electronic screening. Due 
to the large Coulomb interactions, the defects are prob- 
ably incorporated during growth in partially correlated 
positions — oppositely charged donors and acceptors pre- 
fer to sit on nearby sites — and these correlations may 
increase with annealingii^ti^ In Ref. it was shown 
that equilibration of defects during growth or anneal- 
ing leads to an enormous reduction of the typical width 
{{V^) - (V^)2)i/2 of the disorder potential V and to a very 
short correlation length of V, of the order of the lattice 
constant. Ionic screening is thus very effective, whereas 



electronic screening is not. However, the width of the 
disorder potential is still roughly of the same order as 
the Fermi energy so that it cannot be neglected. 

Since the correlation length is so short, a description 
in terms of a delta-function correlated disorder potential 
is reasonable. In this approximation, a scattered carrier 
tends to loose all its momentum information. This allows 
for a relatively simple description of the scattering in the 
semiclassical Boltzmann approach. The spin-exchange 
scattering, though typically weaker than the Coulomb 
scattering, is expected to become important close to the 
Curie temperature Tc, where spin fluctuations are en- 
hanced. A systematic study of the effect of both types of 
scattering on the linear response of DMS and in partic- 
ular on transport would be desirable. For example, the 
resistivity p of (Ga,Mn)As shows a maximum or at least 
a shoulder at r^jgiiLiSiiSiSS whereas the standard Fisher- 
Langer theorjiSi for fluctuation corrections to the resis- 
tivity in ferromagnetic metals predicts an infinite deriva- 
tive of p at Tc- The origin of this weak critical behavior 
is that the resistivity is dominated by scattering events 
with large momentum transfers q ^ 2kp, where kp is the 
Fermi momentum. By contrast, the magnetic suscepti- 
bility x(q) of ferromagnetic metals, of Ornstein-Zernicke 
form)22i2^ diverges only at q = 0. 

As a step towards a comprehensive theory of disor- 
der effects on linear resonse and transport in DMS, we 
present a semiclassical theory for the paramagnetic phase 
of DMS in the metallic regime. Starting from the Zener 
mode&^*^ and semiclassical Boltzmann equations, hy- 
drodynamic equations of motion for the carrier and im- 
purity spin magnetizations are derived in Sec.^ includ- 
ing Coulomb scattering and spin-exchange scattering off 
magnetic impurities. Because of the semiclassical ap- 
proach, these equations hold for small momenta q and 
frequencies lu. The theory is rather general in that it ap- 
plies to both the conduction and the valence band, IITV, 
IT VI, and group-IV host semiconductors, and impurities 
with general spin S. From the equations of motion, the 
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dynamical spin susceptibilities of carriers and impurities 
are derived for small q and lo. The resulting semiclassical 
susceptibility is not of Ornstein-Zernicke form. However, 
this form is presumably restored by quantum effects for 
q of the order oikp- The semiclassical results exhibit the 
detailed dependence on the various sources of scattering. 
We find significant differences between the conduction- 
band (n-type) and valence-band (p-type) cases due to 
the pronounced spin-orbit coupling in the latter. For ex- 
ample, spin diffusion in the valence band is anisotropic. 
On the other hand, we show that semiclassically Berry- 
phase effectsS^iS^ are absent from the linear susceptibility 
even in the valence-band case. 

It would be interesting to study the effect of spin fluc- 
tuations on the electrical conductivity close to Tc in 
DMSj2i This requires the inclusion of quantum effects 
at the scale of kp and thus goes beyond the Boltzmann 
approach. The present theory should be a good starting 
point for this generalization. 

We briefly comment on related work. Sinova et al?^ 
consider the damping of spin waves in the ferromagnetic 
phase in the limit q = within a Green-function ap- 
proach. Disorder scattering is incorporated by assum- 
ing a constant nonzero quasiparticle lifetime. Galitski et 
alm^ derive the local dynamical spin susceptibility close 
to Tc for the strongly localized regime, opposite to the 
case of weak disorder scattering considered here. In the 
strongly localized case the system can be mapped onto a 
disordered ferromagnetic Heisenberg model and GrifRths- 
McCoy singularities are important above Tci2i Qi and 
Zhang^ consider spin diffusion in non-magnetic materi- 
als within the Boltzmann approach. The present work 
goes beyond Ref. ,2^ in that we derive the coupled dy- 
namics of carrier and impurity spins in DMS, consider 
both conduction and valence bands explicitly, and derive 
the dynamical susceptibility. 

As an application we derive the fluctuations of the 
anomalous Hall voltage in the paramagnetic phase in 
Sec, mil In the absence of an external magnetic field the 
average anomalous Hall voltage is zero since the aver- 
age magnetization vanishes. However, fluctuations of the 
magnetization lead to nonzero Hall-voltage noise. Three 
mechanisms of the anomalous Hall effect (AHE) are dis- 
cussed in the literature: skew scattering^^ and side-jump 
scattering^ rely on the imbalance of scattering to the 
right and to the left due to spin-orbit coupling. On the 
other hand, Berry-phase effects"^^ lead to an AHE in the 
presence of spin-orbit coupling even without scattering. 
Since Jungwirth et al^ show that the latter contribu- 
tion can explain the experimental results for DMS in the 
ferromagnetic phase, we also assume this mechanism. 



II. SEMICLASSICAL THEORY 

In this section we present the semiclassical theory for 
the linear response of the carrier and impurity spin mag- 
netizations in DMS in the paramagnetic phase. We first 



derive hydrodynamic equations of motion for these mag- 
netizations and for the carrier magnetization current. 
Some details are given in App. fXl and IbI In App. lUl we 
show that Berry-phase corrections are absent from the 
equations of motion. Then we solve these equations to 
obtain the spin susceptibility. The derivation is carried 
through for both the conduction and the valence band, 
and for arbitrary impurity spin S. We use ?i = A;b = 1- 



A. Hydrodynamic equations, conduction band 

We start with the simpler case of conduction-band elec- 
trons exchange-coupled to impurity spins. Spin-orbit ef- 
fects can be neglected here since the conduction band 
has mainly s-orbital character. This description is appro- 
priate for n-type DMS with nonzero Curie temperature. 
Ferromagnetism in n-type DMS is hard to achieve due to 
the small exchange interaction between electron and im- 
purity spins and is restricted to very low temperatures^ 
We assume a spherically symmetric band ep to avoid 
inessential complications. 

We first briefly motivate the Boltzmann equations for 
the electron density npo-(r), where a = ±1/2 is the spin 
orientation, and for the occupation fraction of impu- 
rity spins with quantum number m of . The Hamilto- 
nian reads 



H = iJkin + J / d'Vm(r) • M(r) 



+ .ge/iB J d3rm(r).Br'+5iA^B J M(r) • B^*, (1) 

where m and M are the electron and impurity spin den- 
sities (oriented oppositely to the magnetizations), re- 
spectively, averaged over microscopic volume elements 
and their coupling is described by the exchange integral 
J = 50 ± 5 meVnm'^.^^ J > (J < 0) corresponds to 
antiferromagnetic (ferromagnetic) coupling. We have in- 
troduced two distinct external magnetic fields BJ;''* and 
B"'^* acting on electron and impurity spins, respectively, 
in order to obtain the linear response of each species sep- 
arately, which will prove useful in Sec. IIIII 

The exchange term is decoupled at the mean-field level. 
We can restrict ourselves to collinear spin configurations 
since the paramagnetic susceptibility is proportional to 
the unit matrix in our spherical model. We choose the 
magnetization direction as the z axis. The mean-field 
Hamiltonian of the electrons and the impurities is then 



Hvin + 9cm I d^rm{r)B, 



Hi = gifiB J d^rM{r)Bi. 



(2) 
(3) 



respectively. In terms of the spin magnetizations fic = 
—gcfJ-B i'm) , /ii = — 5i/iB (M) , the effective fields read 



B„ — B° 



J 



(4) 
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2 A^c- 



(5) 



The single-particle energy of an electron with momentum 
p and spin a = ±1/2 is E^^ = Ep + geMBfiJc- The energy 
of an impurity spin with magnetic quantum number m 
is E]^ = gifiB^nBi. In the absence of scattering, the 
semiclassical equation of motion for the electron density 
riper (r) is given by the Poisson bracket 



= i^rE;, 



(Vpnp,) - (VpE^,) ■ (Vrnp,) 



P"'P(T 



(6) 



with the spin-dependent force — — goMBf Vr-Bc and 
the band velocity Vp = p/wcb, where TOcb is the effective 
mass at the Fermi energy. We use the short-hand nota- 
tion dt for d/dt. With scattering included we obtain the 
Boltzmann equation 



{dt -I- Vp • Vr 



Vp)^ 



f'pa 



5, 



per, 



(7) 



where Sp^ represents collision integrals describing vari- 
ous sources of scattering as discussed below. 

For the impurity spins we define the occupation frac- 
tion of spins with magnetic quantum number m as /,„, 
where J2m f"i ~ ^- corresponding density is 7ii/,„, 
where rij is the density of magnetically active impuri- 
ties. We neglect the contribution of interestitial magnetic 



impurities 
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The Boltzmann equation for the den- 



sity Tiifm is simply dtUifm — Sm, since the impurities are 
assumed to be immobile and purely local. 

We now discuss the collision integrals. The simplest 
one describes disorder scattering of the electrons 



1 



(27r)3 N{0)t 

npcr (1 - "-p'ct)] 



(27r)3 iV(0)7 



<5(ep - ep/ 



^P^)- (8) 



Here, N{0) is the density of states at the Fermi energy for 
one spin component and 1/r is the transport scattering 
rate. Note that there is no change of spin a. 

The next contribution is spin-exchange scattering be- 
tween electron and impurity spins. For this we need the 
transition probabilities between spin states. We write the 
spin operator of the electron (impurity) as s (S). The 
joint spin state is denoted by \am). The matrix elements 
of the exchange coupling are 



2'^<t,i/2^<t' -i/2<5m+i,m' \^S{S + 1) - m{m + 1) 



+ ■^Sa,~l/2^a',l/2^m-l,m'VS{S + 1) - m(m - 1) 



(9) 



Note that only the p' = p contributions to the s^S"^ 
term are taken care of by the mean-field decoupling. For 



p' p this term expresses that carriers can also scat- 
ter off impurities due to the exchange interaction with- 
out flipping the spins. The mean- field approximation 
neglects the discreteness of the impurity spins, which 
this scattering term restores. The transition probabili- 
ties Parn,iT'm' between the states are given by the absolute 
square of the matrix elements, 

^ T'5<T,i/2'^fT' -i/2^m+i,m' [S{S + 1) - m{m + 1)] 



+ ^-f^CT -i/2'5ffM/2'5m-i,m' [S{S + 1) - m(m - 1)] 

(10) 



The collision integral for electron-impurity spin scatter- 
ing can then be written as 



'^P-n J (2^)3 7V(0)r,pi„ 

+ gifJ-BfnBi - €p' - ,gcMBf^'^c - gi^Bm' Bi) (11) 

^ Pam.a' m' [^p'cr' (1 ^pcr) fm' ^pcr (1 ^^p'o"') /m] 

with the spin-exchange scattering rate 1/rspin-'^'* Due to 
conservation of the total spin by the process expressed 
by Eq. (|ll|l . the same collision integral also appears in 
the Boltzmann equation for fm- It is the only scattering 
term we consider for the impurities. 

The scattering processes expressed by iS'^"^ and 5*^?'" 
are not sufficient for a reasonable thermodynamic de- 
scription, however. The reason is that both processes 
conserve the total spin. Thus the homogeneous spin sus- 
ceptibility would be zero. To avoid this problem we al- 
low relaxation of the total spin so that the system can 
approach its thermal equilibrium. This relaxation is im- 
plemented by an additional "spin-flip" scattering term for 
the electrons. Physically, this can be due to the hyperfine 
interaction with nuclear spins^ or electron-electron in- 
teraction in conjunction with spin-orbit coupling in other 
bands. This process is expressed by 



d^p' 1 

(2^)3 iV(0)rfiip 
fip' - gcfJ-Bo'Bc] 



^(ep + goMBcSo 



(12) 



where a = —a. 

The Boltzmann equation for the electrons now reads 

{dt + Vp . Vr + . Vp) np.(r) = S^^^+S'^^+Y, ^^^t 

m 

(13) 



and for the impurities 



(14) 



By summing Eq. H13(l over p, a one easily derives the 
continuity equation dtp + ^r ■ j = for the electron num- 
ber density p = J cPp /{2Tr)^ Upa- and current density 
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j = / d^p/(27r)^ Vp Tipo-. Our main goal is to derive 
corresponding equations for the magnetizations 

= -ffoMB J ^crnp^, (15) 
Mi = -gifJ'Bn.i^mfjn (16) 

m 

and the electron magnetization current 



(17) 



We start with the impurity spins. Multiplying Eq. II14I) 
by m and summing over m we obtain 



7 (27r)3 2^""='P- 



(27r)3 

S'(5+l) Mc , A^(0)T 



3Tspin ScMB 2Tspin ^cMB^i 

iV(0)S'(S'+ 1) 



6t, 



(50MB-B0 - gi^B^i), (18) 



spin 



to linear order in the effective fields and magnetizations, 
cf. App. ^ In the last expression we can identify the 
Pauli susceptibility of free electrons with density of states 
-/V(O) per spin component and the Curie susceptibility of 
non-interacting impurity spins with spin quantum num- 
ber S and density riit^ 



XPauli 



XCu 



A^(0).9eVI 



3T 



(19) 
(20) 



Using these susceptibilities we write 

5(5 + 1) ffi , ^ , 

OtfJ-i = (Mc - XPauli^c) 



3^spin 9t 

1 N{0)T 



2r., 



(Mi - XCurioSi) . (21) 



The rate of change of the impurity magnetization /ij thus 
depends linearly on the deviations of /ie and /Xi from their 
respective equilibrium values, which is quite reasonable. 

Multiplying the Boltzmann equation by a and 
summing over p, a we obtain an equation of motion for 
the electron magnetization, 



5eMB 5eMB 



(22) 



where the force term on the left-hand side vanishes since 
the integrand is a total p gradient. The right-hand side 



can be evaluated similarly to the calculation in App. IXI 
and expressed using xpauii and XCuric, 

• f ^ 5(5+ 1)\ , 

dtHc + Vr • = - ( — 1 — ) (Mo - XPauliSo) 



Tflip 3Ts pm 



1 N{0)T 



2t. 



spin 



— (Mi — XCurieSi) . 



(23) 



To eliminate the magnetization current j^, we derive its 
equation of motion by multiplying Eq. H13() by avp and 
summing over p and a, 



5,^B ' J (2^)3 5 



CTVp (Vp • Vrripo 



d^p 

^ (TVp .goMBC (Vri?c) • VpTlpa 



(27r)3 ^ 
3 .gcMB 4TOcb 

Vr (Mc — XPauli-Bc) 

3 ffeMB 

d^p 

f27r)3 



avJs^: + Sf^ + Y,S^^Z.)- (24) 



The first term — c^tj^/goMB is neglected since it only be- 
comes relevant for frequencies of the order of the largest 
scattering rate. In the second term we have replaced 
WpW^ in the usual way by (5q/5Wp/3, where wf is the Fermi 
velocity. This is valid since np„ has significant r depen- 
dence only close to the Fermi energy. The third term has 
been expanded to linear order in the perturbation and in 
the final step the equilibrium electron density has been 
written as p^°'^ = 2A^(0)mcbw|/3 for a parabolic band. 
Evaluating the integrals, we obtain 



= -D Vr (Mc - XPauli^o) 



(25) 



with the diffusion constant D = w|Ttot/3 and the total 
scattering rate 



1 1 1 5(5+1) 
= - + + -■ 

Ttot T Tflip 4Tspin 



(26) 



Inserting this result into Eq. (|23|l we find the equation of 
motion of the electron spin magnetization. 



9tMc 



spin 



DV^ (Mc - XPauliSo) 



2r. 



1 N{0)T g 

— iMi — XCurioi>ij 



"i ffi 



(27) 



We observe that also the rate of change of Mc is linear 
in the deviations of the hole and impurity magnetization 
from their equilibrium values. The result that 9tMc van- 
ishes in equilibrium must hold in general, not just for a 
parabolic band, as expressed by the Einstein relation. 
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The two equations (|21() and H27|l are coupled both ex- 
phcitly and through the effective fields. They are for- 
mally solved by Fourier transformation in space and time, 



S{S+1) 



Dq- (Mc - XPauli^c 



+ ^™^^(Mi-XCu.cSO, (28) 
S{S +1) gi 

—:; [f^c — XPaulii^oj 

1 N{0)T 



if^i — XCuiic^i) ■ 



(29) 



From these equations we can infer the mean-field Curie 
temperature Tc- In the absence of external fields, 
the static, homogeneous magnetizations satisfy fic = 

-XPauli -/w/SeSiMB and fii = -XCuTieJlJ-e/9e9ifJ'B "^'^^^ 

nonzero solutions at the Curie temperature^i^iSL^SiSS 

T,= ^^^±2lNiO)j'ni. (30) 



B. Hydrodynamic equations, valence band 

We now derive hydrodynamic equations for valence- 
band holes exchange-coupled to impurity spins, relevant 
for p-type DMS. The case of spin quantum number S = 
5/2 corresponds to substitutional Mn in GaAs. The main 
complication here is the presence of spin-orbit coupling. 
We employ a 4-band Kohn-Luttinger Hamiltonian^i^Si^ 
in the spherical approximation, which is the simplest one 
incorporating the relevant physics. In the absence of 
magnetic impurities the Hamiltonian reads^M^ 

[ill + 572/2) P - 272 (k • j)2] (31) 

with Kohn-Luttinger parameters 71, 72 and the angular 
momentum operator j of the holes, which in this subspace 
can be written as a 4 x 4 matrix and has the Casimir op- 
erator j • j = 3/2(3/2 -I- 1). Since the split-off band is ne- 
glected, this description only applies to semiconductors 
with sufficiently strong spin-orbit coupling. The eigen- 
states of -ff at k are characterized by the quantum num- 
ber j = ±1/2, ±3/2 of k • j, where k is the unit vector in 
the direction of k, i.e., the spin quantization direction is 
k. We restrict ourselves to the heavy-hole band, which 
is justified for energies close to the band edge because of 
the much smaller density of states of the light-hole band. 

We introduce the eigenstates | j)k of k • j with eigen- 
values j. We denote the spin eigenstates with respect to 
a fixed quantization axis z by The former can be 
expressed in terms of the latter by means of a rotation 



m spm space 



43.44 



|J), 



(32) 



where j^, are spin operators and 9 and are the polar 
angles of k. 



The states |j)k can be expressed in terms of eigenstates 
of orbital angular momentum 1 (with 1-1 = 2) and spin s 
with the help of Clebsch-Gordon coefficients. One then 
easily finds that all 4 x 4 matrix elements of equal 
the corresponding matrix elements of j^/3. The same 
holds for the x and y components because of symmetry 
so that s = j/3 holds as an operator identity in the heavy- 
hole/light-hole subspace.'^"' Consequently the heavy-hole 
states {j — ±3/2) are eigenstates to k-s with eigenvalues 
±1/2. However, the heavy holes alone do not form a spin 
doublet since the matrix elements of = Z^/3 all vanish 
in the two-dimensional heavy-hole subspace — single spin 
flips cannot change the total angular momentum from 
+3/2 to —3/2 or vice versa. 

The band energy of the heavy holes is ep — (71 — 
272) p^/2m. Together with the Zeeman energy their total 
energy is 



Ep^ = Cp + .ghMB ^ cos 6 Bh, 



(33) 



where i?h is the effective magnetic field and 9 is the polar 
angle of p with respect to the field direction z. Without 
scattering the equation of motion for the hole density 
reads 

dtupj = -{npj,E^^} 

j p 

= gtifJ-B 7; cos 9 VrSh • VpWpj • Vr^pj 

6 mhh 



.giiMB I sin 6*511 - • VrTT-pj, 

3 p 



(34) 



where rrihh = '^1/(71 — 272) is the heavy- hole effective 
mass. This suggests to define the velocity as 



.ghA^B o sm6'Bh -. 

mhh 3 p 



(35) 



Note that the second term is explicitly of first order. We 
should use this velocity in the semiclassical equations. 
However, we find the contribution from the second term 
to vanish to first order. The reason is essentially that 
we have to evaluate all other factors in equilibrium due 
to the explicit J5h. This result is proved together with 
the absence of Berry-phase contributions in App. [CI We 
thus drop the second term in Eq. (|35|l . 

We now turn to the derivation of the Boltzmann equa- 
tion for the holes. Analogously to the conduction-band 
case we have 



{dt + Vp . Vr + Fp, . Vp) np, (r) = 5^;^ + Y.^\ 



ispm 
pjm 



(36) 

with the force Fpj = — .ghMB (j/3) cos6'Vri3h for the 
holes and 



dPp 
(2^ 



Et^spii 



(37) 



6 



for the impurities. The disorder scattering in- 
tegral contains the matrix elements k(jb')k' — 



{j\e 



The spin operators are 



4x4 matrices in the projected subspace. For heavy holes, 
explicit evaluation gives the transition probabilities 



c(j|/> 



cos 



sm 



a 


sin^ 


a 






2 




2 


a 


cos^ 


a 






2 




2 



(38) 



where — ±3/2. Here, a is the angle between the 
vectors k and k'. The collision integral for disorder scat- 
tering of heavy holes reads 



P^' / (27r)3 ^ 7V(0)t 



<5 f Cp + ffhAiB I cos 8 Bh 



-V -5hAtB — cos6''Bhj |pO|j')p'| (n-p'j' - ripj). 

(39) 

Note that for forward scattering (a 0) we get predom- 
inantly j' — j, whereas for backscattering (a ^ tt) we 
find predominantly j' = —j. 

Due to the k-dependent quantization axis the quantum 
number j is not conserved even by pure disorder scatter- 
ing due to the Elliott- Yafet mechanismi^ This scattering 
takes a hole of momentum k and quantum number j into 
a state of momentum k' under conservation of spin. How- 
ever, its spin state is no longer an eigenstate at k'. In the 
semiclassical approximation it assumes possible magnetic 
quantum numbers j' with probabilities |pO|j')p'P- 

For the hole-impurity spin scattering we need matrix 
elements of s • S. The transition probabilities are 



= \p{jm\s ■ S|j'm') 



P3in,p'fm' 

^1(1 Ip01j + I/)p'I' '5™+i.„. [S{S + 1) - m(m + 1)] 
+ I \p{j\r\j')p'fS^^i,m' [S{S + 1) - m(m - 1)] 

+ \p{j\f\j')p'\^Smrr,m^^ (40) 

and the resulting collision integral reads 



gifiBmBi - Ep, - ffhAiB-^ cos 9' By, - gi^Bm'Bi 



X R 



pjm,p'j'm' 



(41) 



Since the two collision integrals already include spin re- 
laxation we do not introduce an additional spin-flip term. 



We now derive hydrodynamic equations for the hole 
and impurity spin magnetizations. Some details of the 
calculations are shown in App.^ The hole and impurity 
spin magnetizations are 



(2^) 

J 

Mi = -giA'B^-i "T- /m 

and the hole magnetization current is 



/(gaEi-^-p.^ (42) 



(43) 



-5hMB 



^ I coseiVpTipj. (44) 



(27r)3 ^ 3 



We start with the impurity spins. In analogy to the 
conduction-band case we obtain 



S{S + l) ( I „ 



5 iv(o)r 



{pi - XCurio^i) , (45) 



where we have again identified the Pauli susceptibility 
XPauii = ^(0)5hA'B/2 and the Curie susceptibility. This 
result is of the same linear form as for the conduction- 
band case. Note, however, the presence of a factor of 1/3 
multiplying the Pauli susceptibility, which is absent for 
the conduction band. This factor is easily understood by 
calculating the static, homogeneous spin susceptibility of 
heavy holes in the absence of impurities. For the static 
susceptibility we can assume the holes to be in thermal 
equilibrium. 



'pj - nF ( Cp - M + .9hMB I cos 6 Bh 



= npiep - m) + (fp - M)ffhMB ^ cosOB^, (46) 

where n^p {E) = npiE) [npiE) — 1] is the derivative of 
the Fermi function. To linear order in B^ we thus find 



= -ffh^B 



(27r)3 ^ V3 

3 



2..2 I ^(Zj cosHBhu'^piep-^) 



1 N{0)glf^l _ I 

2 2 -^h — gXPauli-Dh- 



(47) 



The extra factor stems from the angular integral over 
cos^ 9 and is thus due to the heavy holes not being simple 
spin-1/2 fermions, as discussed above. 

We also obtain the equation of motion for the hole 
magnetization. 



7 



——J 2_, 3 cos6'5hAiB g cose* (VrSh) • VpHpj 

j 



ffhMB ffhMB y (27r) 



,,,3 ...B y(2.)3E^O«H"^^"?"^"j' ^^^^ 

to first order. Similarly to the calculation in App. ^ we obtain for the right-hand side 

To eliminate the magnetization current we consider its equation of motion. The left-hand side is 



(2^ 3 cosSvp (vp • VyUpj) - j 2l 3 cos^VpghMB ^ cos6' (Vr-Bh) • Vp^p^. (50) 



I 



The first term is again neglected. In the second we have 
to be more careful because of the explicit angle depen- 
dence. For the conduction band the factor Wp /3 is ob- 
tained by assuming np^ to be the equilibrium distribu- 
tion in a constant Zeeman field. The integral over the 
direction of p is then easily performed. Since we obtain 
a term linear in Vr/^e, corrections would be of higher 
order. For the valence band we also assume a constant 
Zeeman field, leading to Upj = Up^ + {j / 3) cos 6 An{p) . 
Thus the second term in Eq. (|50|l becomes 



d^p 
(2^ 



^ cos Op (p ■ Vr^ COS 6 An{p) 



1/15 
X I 1/15 I VrAn[p(0]. (51) 
1/5 

In the same approximation we find —fJ-h/ghfJ-B 
^(0) / Ej 073)^An[p(f)]/3 so that this term is 



1/5 , 
^2 I 1/5 Vr-^ 

3/5 / 5hMB 



(52) 



The third term in Eq. I|50(l is straightforward to evaluate 
to first order, 



N{0)gi,m ..2 
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1/5 
F , 1/5 I VrSh, (53) 
3/5 



for a parabolic band. Again, the result holds in general 
due to the Einstein relation. Altogether the equation of 



motion for the magnetization current is 



1/5 

Vp I 1/5 

3/5 

d^p 



1 



■ XPauli Bh 



(27r)= 



E 



|cos0vp(5^^ 



E<^spii 
Pi" 



(54) 



Evaluating the integrals we finally obtain 



3/5 
3/5 
9/5 



— g XPauli^h 



(55) 

where we have introduced the diffusion constant D — 
VpTtot/S with the total relaxation rate 1/rtot — 1/(2t) -t- 
5S{S + l)/(72Tspin). The spin diffusion in the valence 
band is thus anisotropic. Compared with the result H25|) 
for the conduction band, diffusion along the direction of 
the effective field is enhanced and diffusion in the trans- 
verse directions is suppressed. The origin of this inter- 
esting effect again lies in the momentum dependence of 
the quantization axis due to spin-orbit coupling in con- 
junction with the projection onto heavy holes: Consider, 
for example, heavy holes traveling exactly along the x di- 
rection. In the Hilbert subspace of these holes all matrix 
elements of and vanish so that these holes cannot 
carry any spin magnetization pointing in the z direction. 
For holes with momentum p pointing mostly but not 
fully in a transverse direction the contribution to spin 
transport is still suppressed. 

Inserting our result for the current into Eq. H49(l we 
obtain the equation of motion for the hole magnetization. 



1 75(5 + 1) 



5t 



ISOTspin 
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- D 



5 dx'^ 



5 dy^ 



36r, 



1 N{0)T 5h , 



9 ^ 

5 dz^ 



XCuricSi) 



(56) 



This equation is of the same general form as for the con- 
duction band, the main differences being the reduced 
Pauh susceptibiUty and the anisotropic spin diffusion. 

As in the conduction-band case, Eqs. and are 
coupled by the effective fields Bh = B^^^ — ( J/ghSiMl) Mi 
and Bi — Bf^^ — ( J/ghffiA^B) A^h- Fourier transformation 
yields 



1 75(5 +f) 3ql + 'Sql + 9q^_ 

1 ^ + D 

5r fSOTspin 5 



X ( Mil ~ 3 XPauli-Sh 

— (Mi - XCuric-Bi) , (57) 



36Tspin rii gi 

fj i ( 
— Mh 
5h V 

5 iv(o)r 



—rZ Mh — -r XPauli£*h 



36t, 



(Mi - XCuricSi) ■ (58) 



spin "-1 



For uo — Q, q = we find finite solutions at 



T = Tc = 



5(5 -fl) 
fS 



N{0) Hi. 



(59) 



The Curie temperature of holes is reduced by an extra fac- 
tor of f/3 compared to the conduction-band case. This 
factor stems from the same factor in the Pauli suscep- 
tibility, which is ultimately due to spin-orbit coupling, 
as discussed above. On the other hand, for typical host 
materials the density of states is much higher for the 
heavy holes than for conduction-band electrons and the 
exchange integral J is also much larger, enhancing Tc in 
p-type materials. 

We have so far ignored the possible effect of Berry- 
phase contributions. As noted in SecO] these have been 
found to be important for the magnetic response of p- 
type DMS, where they contribute to the anomalous Hall 
effect in the ferromagnetic phase.-"* ft is thus necessary 
to check whether they contribute to the hydrodynamic 
equations. In App. [21 we show that they do no con- 
tribute to linear order in the effective field. Berry-phase 
contributions are expected in higher orders, though. 



C. Susceptibilities 

With the help of the hydrodynamic equations we now 
derive the linear response of the carrier spin and impurity 
spin magnetizations to external fields coupled to these 
magnetizations. It is useful to solve the general problem 



of the Fourier-transformed hydrodynamic equations 

-iujflc = (Mc - "XPauliSc) 

+ i?ci - (Mi - XCuricBi), (60) 



9i 



-iuJfJi = Ric — (Mc - "XPauli-Bc 

9c 



iv(o)r 

-Kii (^i - XCurioiSi ) , (61) 

Hi 



where 



B, = B^„ 



2" Mi 

5c.giMB 



9 Me 5 

5c5iMB 



(62) 
(63) 



which contain the special cases of the conduction band, 
Eqs. (|2H1) and (j^Hl), and the valence band, Eqs. ((23 and 
(|58|l . For the valence band, the subscript "e" should 
of course be replaced by "h". Solving this system of 
equations we find 



Mi/5iMBy detM Ai, ) \g,fi^Bl^^ 



with the determinant of the coefficient matrix 



(64) 



det M = Qni{ibjf + iuj [2aN{Q)JniRic - GriiRq 

+ 2N{0)JniRciS{S + 1) - 6A^(0)Ti?ii] (65) 
4- iV(0) {RciRic - RiiRq) [aN{0)J^niS{S + 1) - 6T] 



and 



Aec = -3a[ic^nii?q + 7V(0)T(i?eii?ic-i?iii?q)], (66) 
^oi = an; [3zcji?ie + iV(0) J (i?eii?ic - i?iii?q) 

X 5(5+1)], (67) 
Ai, = ni[2iujRei + aN{0)J{R,iRi, - RiiRq)] 

X 5(5+1), (68) 

Ai = -2ni[i?eii?ie + i?n(jc^-i?q)] 5(5 + 1). (69) 



The magnetization becomes singular at 

T^T,= -SiS + l)N{0) J^nu 
6 



(70) 



in agreement with our earlier results.^® We now assume 
u and T — Tc to be small compared to the rates i?q, Rd, 
-Rio, -Rii but do not make any assumption about T ~ 
vs. w. Then we find 



Xcc Xci 

Xio Xii ) \ B-?""^ 



(71) 



with the susceptibility matrix 
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Xio Xii 



2iV(0) 5(5 + 1) (i?eii?ie - i?iii?q) A^l 



X lUJ 



[6i?q - 3aiV(0) Ji?ie - 2N{0)JR^iS{S + 1) + aN^{0)J^S{S + l)Rn] 



aN'^{Q)J^S{S + 1) {RciRic - RiiRq) 



r 



^ I 2 ( aN{Q)J 
9c 1 2 



V 



a7V(0)J 



(72) 



^ aN{0)J \ 

gf 



Since we cannot apply different fields to tfie carrier spins 
and the impurity spins, the physical susceptibility of the 
carrier spins is Xcc + Xci, while the susceptibility of the 
impurity spins is Xic + Xn- The total susceptibility de- 
scribing the response of the total magnetization is 



Xtot = Xcc + Xci + Xio + Xii- 



(73) 



Note that this physical susceptibility is always para- 
magnetic since the components of the matrix factor in 
Eq. lO combine to {gcaN{0)J/2 - g^Y . In the static 
case w = all four components are of Curie form. We al- 
ready see that the dimensionless parameter —aN{0)J/2 
has a special meaning: It is the ratio between the av- 
erage electron spin and the average impurity spin in an 
applied field, regardless of whether the field acts only on 
the electrons, on the impurities, or on both. 

We now consider the special case of conduction- 
band electrons. Inserting the appropriate factors from 
Eqs. (|28|l and H29(l . we obtain the susceptibility matrix 



X = iV(0) 5(5 + 1)m| 



6r,pi„ + 25(5+1) 
25(5+1) 



(1 - iV(0) J/2)2 



2/Tflip + Dq^ 



V 2 

/ 2/iV(0)J^' 



9c 



V 



V 2 

N{0)J 
^9c9i — — 



'9c9i 



Tc 

N{0)J \ 



9i 



(74) 



with = [5(5 + l)/6]iV(0) J^n;. This susceptibility 
describes the linear response of an n-type DMS. The same 
result would be obtained for a simple model of spin-1/2 
holes, which is sometimes employed in the literature. 

Note that the only q dependence appears in the co- 
efficient of Lu. This is quite different from the standard 
Ornstein-Zernicke formSSiS^ of the susceptibility. We dis- 
cuss this point further below. The only typical length 
scale in x is = \/ -DTflip/2. This is the relaxation length 
of the total spin since the total spin relaxes with the spin- 
flip rate l/rflip. In the semiclassical approximation 
does not show any critical behavior at T^,. 

We show below that N{0)J <C 1 for the valence band. 
For the conduction band the density of states A^(0) is 



smaller than for the valence band and the exchange inte- 
gral J, which for the conduction band is predominantly 
due to onsite Coulomb exchange, is also smaller so that 
N{0)J would be very small for n-type DMS. 

The susceptibility also describes the magnetic excita- 
tions. Their dispersion is obtained by equating the de- 
nominator to zero and solving for oj. We see that these 
modes are diffusive with relaxation rates 



25(5+1) 



X — iuj — 



f N{0)J 



6ts, 



25(5+1) 



(1 - jV(0) J/2)^ 
2/Tflip + Dq^ 



T-T, 



(75) 

The rate A is always positive for T > Tc, as required for 
exponentially decaying excitations, and is smallest for 
q = 0. The q dependence is controlled by the total-spin 
relaxation length ^c- In the semiclassical approximation 
A goes to zero for T — > for all q simultaneously, but 
see the discussion below. 

We now consider the case of valence-band holes. In- 
serting the appropriate parameter values from Eqs. H57|l 
and H58|l we obtain 



18 



5(5 +1)7V(0)m| 



5(5+1) (1 - 5A^(0) J/6)^ 
^ 15 

fN{0)J 



Rci 



( „2 ^^(0)^ 



V 



7V(0)J 
-5o.gi — 7. — 

6 



'9c9i 



Tc 

N{0)J \ 



6 



9i 



(76) 



with the Curie temperature Tc = [5(5+ 1)/18] N{Q) J^n; 
and 



■ — RhiRih — RiiRq 



36re 



spin 
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36r„, 



5Tt, 



+ D 



Ml + H 




(77) 



This susceptibility applies to p-type DMS. The result 
is of nearly the same form as for the conduction band. 
The only differences except for simple rescaling is that 
not the parameter — iV(0) J/6 itself but — 5iV(0) J/6 ap- 
pears in the q-dependent term and that the diffusion is 
anisotropic. 

For both the conduction band and the valence band the 
susceptibilities depend on q only through the coefhcient 
of the frequency u>. The static susceptibility {u> = 0) is 
thus independent of q in our approximation. This would 
mean that the instability appears simultaneously at all q. 
The tendency of the system to become ferromagnetic is 
not found within the semiclassical Boltzmann approach 
since it does not incorporate physics at large momenta 
q ^ kp. We expect the most important effect for q kp 
to be the q dependence of the Pauli susceptibility.^^ In- 
serting this dependence by hand, we obtain an additional 
term of the order of +q'^/kp in the denominator, which 
makes the instability first appear at q = 0, leading to fcr- 
romagnetism. A rigorous evaluation of the susceptibility 
at all momenta requires a fully quantum-mechanical cal- 
culation, e.g., using the quantum Boltzmann equation. 
We leave this as work for the future. 

One could think that a ferromagnetic interaction be- 
tween the carriers themselves introduces a new length 
scale and might therefore introduce a q^ term into the 
denominator of x- In our approach such a ferromagnetic 
coupling between the carriers, say holes, leads to an ad- 
ditional term in the effective field, 



By, 



Bh 



J 



.9h.giMB 



K 
5hM| 



(78) 



with K > 0. The derivation can be carried through. 
We only show the resulting susceptibility for the valence 
band (the conduction-band result is analogous), 



X 
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S{S + l)NiO)^il - iu 



S{S+1) (1 - 5iV(0) J«/6)^ 
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N{0)J, 
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(79) 
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with K = N{0)K/6, J« = J/(l - k), andSi^a 

5(5 + 1) N{Q)J^ni 
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1 - K 



(80) 



The Curie temperature is enhanced by the Stoner factor 
(1 — k)^^. _Rq is still given by Eq. H77|l . If the ferromag- 
netic interaction becomes so large that N{0)KlQ = 1 



then the hole system orders ferromagnetically at T = 
even in the absence of any impurity spins (Stoner 
instabilitji2&) and our approach breaks down. However, 
for DMS K is small. ^ The same result is obtained by in- 
troducing an appropriate Landau parameter Fq — —k 
into Fermi liquid theoryiiiSr We see that the inclusion 
of a carrier-carrier ferromagnetic exchange interaction 
changes the susceptibility quantitatively but not its func- 
tional form. In particular, it does not introduce a 
Ornstein-Zernicke-type q^ term. 

Let us estimate the parameter N{0)J/6: For a para- 
bolic band 



7V(0) 



27r2 



27r2 



(STT^n) 



1/3 



(81) 



where n is the carrier density. For Gai-ajMn^jAs with 
X = 0.05 and p = 0.3 holes per manganese atom we get 
iV(0) « 7.26X lO-^meV"^ nm-^. On the other handjii^^ 



J « (50 ± 5) meV nm"* so that 

iV(0) J 
6 



0.0061. 



(82) 



The parameter is thus small. However, we emphasize 
that the derivation is valid for general N{0)J. The small 
value explains why the hole contribution to the magne- 
tization is small compared to the manganese one^ii^ 
Also, for antiferromagnetic coupling J is positive so that 
the hole and impurity spin magnetizations are opposite 
in sign, in agreement with experiments^^iiLiSiia 

The above estimate of A^(0) relies on the spherical ap- 
proximation and on the omission of the light-hole band, 
which are not well justified at the hole concentration used 
here. A realistic Slater-Koster tight-binding description 
of the unperturbed valence band^ gives a density of 
states per spin direction of 1.18 x 10~^ eV~^A ^. The 
dimensionless parameter A^(0) J/6 « 0.0099 is thus some- 
what increased by assuming a realistic band structure. 

Equation (|77|) shows that the typical length scale of 
X is = v^5£'Ttot/2, which corresponds to in the 
conduction-band case. The time appearing in should 
thus be the relaxation time of the total magnetization. 

The magnetic excitations are again diffusive modes. 
Their relaxation rates are 



X = iuj = 



f 5(5+1) (mi 



5(5+1) (1 - 5iV(0) J/6)^ 



T-T, 



15 



Ra 



(83) 

qualitatively similar to the conduction-band case. 

We propose to measure the magnetic susceptibility 
in the paramagnetic phase at small q and lu for vari- 
ous DMS. This should allow to test the non-standard 
functional form of our result. In particular, such an ex- 
periment should look for the anisotropic spin diffusion 
in p-type DMS. Studying samples with similar concen- 
trations of magnetic impurities but different concentra- 
tions of nonmagnetic scatterers introduced by codoping— 
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would allow to change the scattering rate 1 /r while hold- 
ing l/rgpin and the mean-field nearly fixed. 



III. ANOMALOUS HALL- VOLTAGE NOISE 



and current density j = e / d'^p/ (Stt)^ J^j '^pj '^pi- From 
a similar evaluation as in Sec. Ill Bl we recover the Drude 
conductivity to order zero in the impurity spins S. 

The first contribution to the Hall current is found at 
first order in S 



In this section we apply the semiclassical theory to the 
derivation of the voltage noise in the transverse direction 
in the paramagnetic phase. The average anomalous Hall 
voltage vanishes for T > due to the vanishing average 
magnetization. However, fluctuations in the magnetiza- 
tion are present and are in fact critically enhanced as Tc is 
approached. This leads to fluctuations in the anomalous 
Hall voltage, which we derive in the following. Following 
Ref. I24I we consider the Berry-phase contribution to the 
anomalous Hall effect for a p-type DMS. 

The fluctuations in the Hall voltage are governed by 
the correlation function of the effective magnetic field 
acting on the hole spins. This correlation function is 
closely related to the impurity-impurity spin susceptibil- 
ity Xii evaluated above. Typical Hall-bar samples are 
much larger than the spin relaxation length ^h- Hence, 
we can restrict ourselves to the homogeneous component, 
q = 0. Fluctuations with nonzero q cancel out in the 
macroscopic voltage measurement. On the other hand, 
the frequency dependence of the xn is important since lu 
can become larger than T — Tc close to the transition. 

We describe a p-type DMS in the metallic regime by 
the Hamiltonian 



H = i/kin + ^'^iS • S — eE • r, 



(84) 



where s is the hole spin operator, S is the averaged (q = 
0) impurity spin and E is a homogeneous, static external 
electric field. The external magnetic field vanishes. The 
kinetic Hamiltonian is i?kin = {I /2m) [(71 + 572/2) — 
272 (p-j)'], as in Eq. IHJ. 

The anomalous Hall conductivity has been derived by 
Jungwirth et aZ.^^ The derivation is similar to the one 
in App. [ni The exchange and electric-field terms are 
treated as small perturbations. The equation of motion 
of r, Eq. (jCip . can be rewritten as24 



r = VpEpj - eE X n 



2lm{Vpu\dtu) (85) 



with fl = Im (Vpu| X |Vpu) and the heavy-hole energy 

Epj = 1^(71 - 272) + ^ j p • S - eE • r, (86) 
2m 3 

up to first order in E and S. 

The charge response is derived from the Boltzmann 
equation (9t-|-r- Vr-l-p- Vp)npj = SpJ. We restrict our- 
selves to nonmagnetic disorder scattering, assuming the 
disorder scattering rate 1 /t to be large compared to the 
spin scattering rate l/Tgpin, since inclusion of the latter 
would only complicate the notation without introducing 
new physics. For the anomalous Hall effect we are con- 
cerned with the charge density p = e J (fip/ {^tt)^ ripj 



.^1-, e Jrii ^ 



1 



71 - 272 372 



where kp is the Fermi wave number in the heavy-hole 
band. A homogeneous charge distribution has been 
assumed to obtain this result. In the limit of large 
heavy- hole/light-hole mass ratio Whh/wih S> 1 the Kohn- 
Luttinger parameters satisfy 71 — 272 ^ 72 and we obtain 
the simpler result j^^^ — e^mhY^Jn-i/ {Air'^kF) E x S. Note 
that the first-order contribution is purely transverse — 
there is no anomalous contribution to the longitudinal 
resistivity to this order. The anomalous Hall current 
density is, to first order, jah = f ah E x S with the unit 
vector S in the impurity-spin direction and^ 



47r2fcp 



In the paramagnetic phase the average spin magneti- 
zation vanishes and, therefore, the average anomalous 
Hall current also vanishes. However, its ffuctuations 
(jAH • Jah) do not. We write Jah = ctE x S, where 
— ctah/S'. Thus 

(jAH-jAH) =^'((ExS) (ExS)) 

= d^(E^{S-S)-Y,E^Ep{S"SP)\. (89) 

^ a/3 ^ 

In the paramagnetic phase this gives 

(jAH(i) •jAH(O)) - 2~a^E^ {S\t)S\0)). (90) 

The time-dependent correlation function can be ex- 
pressed by the impurity-impurity part Xii of the suscepti- 
bility in the p-type case with the help of the ffuctuation- 
dissipation theoremj^fe 



dte-*"*(S'^(r,t) S'^(0,0)) 
^ 2rimxr(r,c^) 

9 9 9 : 



2 Im xfj^(r, ijj) 
gff^lnf (1 - e'-n 

(91) 



where the final expression is valid for lo <ti T. 

We now evaluate the correlation function of the anoma- 
lous Hall voltage C/ah between the front and back sides 
of the relevant Hall-bar region shown in Fig. ^ Since 
Coulomb interaction suppresses charge fiuctuations the 
current density is assumed to be homogeneous; devia- 
tions are only expected to occur at frequencies of the 
order of the plasma frequency. We can then write the 
anomalous Hall voltage as C/ah = Lyjj^^/ajj, where 
ap) = e^rii^T /ra\i\^ is the Drude conductivity. If the elec- 
tric field is applied in the x direction, the anomalous Hall 
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The ratio of conductivities is 



FIG. 1: Geometry of the relevant section of the Hall bar. 



current density in the y direction is — —BES^ . Due 
to homogeneity we can average the current over the sam- 
ple volume. The correlation function of the voltage is 
then 



(C/AH(i)C/AH(0)) 
^2 1 



/ d^rd^r'(5"(r,i)5"(r',0)) 



n-2 r 2 r 2 

= 4 I (^'(r, t) S'iO, 0)). (92) 

Taking the Fourier transform, expressing the electric field 
by the voltage applied to the relevant sample region, E — 
U/Lx, and inserting Eq. (|91|l we obtain 

(UahUah)^ = J dte-''^\UAHit)UAHm 
_ Ly ^2 2Imxr(q = 0,w) 



gftJ-lnf (1 



(93) 



Assuming 1/r > l/rspm and N{0)J < 1, Eq. gives, 
to leading order in u, 



3N{0) gffxl T,pi, 



55(5+1) 



/7V(0)jy fT-T, 



) m 



For cj <C T this leads to 
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spin 



(95) 



so that the noise spectrum is independent of uj for small 
Lo. Close to the Curie temperature the integrated noise, 
(^ah) = (t^AHC/AH)wAti; with the detector bandwidth 
Ao; = 27rA/, satisfies 



^AH/ ^ 

t72 al LIL^ 
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(97) 



with the Fermi energy Ep = A:|,/2mhh- The factor n{/n\^ 
lies in the range 1 ... 10, the ubiquitous factor iV(0) J/6 
drops out of the final result, and 1/EpT has to be rea- 
sonably small for our metallic picture to apply. The final 
dimensionless expression for the integrated noise is 



ah) 



27 



C/2 



1 



1 



5 \nh EpT J L^Lz n-i \T — Tc 



'spin- 

(98) 

This contribution to the noise is critically enhanced as 
the Curie temperature is approached. In a homogeneous 
system it should diverge at Tc but real DMS are, by their 
very nature, disordered and the transition is broadened 
by macroscopic inhomogeneity of Tc- Furthermore, the 
effect strongly depends on the length of the relevant 
region of the Hall bar in the electric-field direction, being 
large for small L^- It is more weakly enhanced by a 
small sample thickness and by a large sample width 
Ly across which the voltage is measured. The effect is 
also increased by strong compensation (nh <C rii) and in 
samples showing bad metallic behavior (small Epr). 

The anomalous Hall-voltage noise is in competition 
with the thermal (Johnson-Nyquist) voltage noise^S 
which in integrated form is {U'^-^^ = 2TR/S.ijj/-k = 
2T{Lyl (TdL^L^^lo j-K . The two contributions can be 
experimentally distinguished by their different temper- 
ature and voltage dependences. The anomalous Hall- 
voltage noise (f/An) proportional to the applied volt- 
age squared, whereas the thermal voltage noise is inde- 
pendent of voltage. 

Besides being an interesting physical effect, measure- 
ment of the anomalous Hall-voltage noise would provide 
an independent approach to the impurity-spin suscepti- 
bility and to important experimental parameters, such as 
the compensation fraction n^^/n^ with respect to the den- 
sity of magnetically active impurities. The Hall-voltage 
noise would also provide a new way to determine the 
Curie temperature. More generally, such experiments 
would test the applicability of the semiclassical theory to 
DMS. 2^ It may also be interesting to study the anomalous 
Hall- voltage noise in conventional itinerant ferromagnets 
such as iron. 



IV. CONCLUSIONS 

A semiclassical approach based on Boltzmann equa- 
tions for electrons or holes and impurity spins has been 
used to derive hydrodynamic equations of motion and 
spin susceptibilities of diluted magnetic semiconductors 
(DMS) in the paramagnetic phase. This theory gives the 
leading frequency and wave-vector dependence at small 
Lo and q. Our results apply to p-type and n-type DMS, 
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to III-V, II- VI, and group- IV host semiconductors, arbi- 
trary impurity spin quantum number S, and ferromag- 
netic or antiferromagnetic exchange couphng J of carrier 
and impurity spins. While the form of the equations of 
motion is easy to understand, the susceptibihty has a 
nonstandard q dependence, which only appears in the 
frequency-dependent term. Thus the semiclassical diffu- 
sive dynamics does not lead to any q dependence of the 
static susceptibility. Such terms are expected to be in- 
troduced by physics at the much larger momentum scale 
of the Fermi momentum kp . 

Spin-orbit coupling in the valence band leads to quali- 
tative differences in the susceptibility of holes compared 
to electrons. The first difference is a suppression of 
the mean-field Curie temperature of p-type DMS com- 
pared to n-type DMS by a factor of 1/3, which can be 
traced back directly to the momentum dependence of the 
spin quantization axis in the presence of spin-orbit cou- 
pling. On the other hand, the Curie temperature in p- 
type DMS is enhanced by the typically larger density of 
states and exchange coupling. The second difference is 
the anisotropic spin diffusion in the valence band, which 
is apparent in the equation of motion of the hole magne- 
tization and also makes the q dependence of the suscep- 
tibilities anisotropic. The anisotropic diffusion is due to 
the fact that holes moving in a direction perpendicular 
to the magnetization or effective field have vanishing ex- 
pectation value of the spin in the magnetization direction 
and thus do not contribute to its transport. 

The results have been applied to evaluate the noise in 
the anomalous Hall voltage in DMS, which is governed 
by the impurity-spin susceptibility at small frequencies 
and momentum q — s- 0. Unlike the average anomalous 
Hall voltage this quantity does not vanish in the param- 
agnetic phase and is even critically enhanced close to Tc- 
The noise gives an independent experimental approach 
to the impurity-spin susceptibility. We have derived the 
detailed dependence of the signal on the impurity and 
hole concentrations and on the sample geometry. 
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APPENDIX A: HYDRODYNAMIC EQUATIONS, 
CONDUCTION BAND 

In this appendix we collect a number of calculations 
pertaining to the conduction-band case. The derivation 
of the hydrodynamic equations in Sec. Ill Al requires the 
evaluation of various integrals over the collision terms 
Sp^, 5pP'^, and 5pJP. We do not show all evaluations 
but only present a few to clarify the method and approx- 
imations used here. 

The first integral we need is 

^ ^ crm 

which appears in the equation of motion (|18|l of the impu- 
rity spin magnetization. We divide the collision integral 
into three terms, 

cspin cspin,0 _i_ cspin, + l i cspin, — 1 ( 

corresponding to m' = m (no spin flip), to' = m -I- 1, and 
to' = m — 1, respectively. The first contribution is 

J (27r)3 P^™ 4Ar(0)r3pi„ 

^ ' am 

as can be seen by renaming p <-> p' in the term with Tip^. 
The other two contributions can be treated together as 



X m [S'(5' + 1) - m{m ± 1)] [yip'zp (1 - np±) fm±i - np± (1 - 7ip'=p) /„] . (A4) 



We write = 1/(25 + 1) + A/„„ where Y.m ^fm = 0, in A/„, 
and divide the integral into terms of zero and first order 

^ ' am 
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In the zero-order term we expand the delta function in 
Bg, Bi, and write all terms strictly in first order. This 
allows to perform the integrals, 



.(0) ^ 



4t ■ ^ 



2STT 



(27r)3 



'P'T 



(2^ 



>g(^+l) Me 
6rspin gcMB 

^ A^(0)5'(5 + l) 



(gcA^BSo - gi^iBBi). 



(A6) 



We have used partial integration in the last term. The 
term S*^^) in explicitly of first order in Afm so that all 
other factors are to be evaluated in field-free equilibrium, 



7V(o)r 



^ ' spin 



t[^(5+l)-m(m±l)](A/„,±i-A/„ 



(A7) 

In the sum we replace m by m =p 1 in the term containing 
fm±i- If we still sum over m from —S to 5*, we expect 
additional contributions at both ends, but these vanish 
due to the factor S{S + 1) — 171(771 ± 1). Thus we obtain 
Ed) = JV(0)T/4r,pi„ j:j-m ± 3m^ T S{S + 1)] A/,„. 
This expression obviously simplifies when the contribu- 
tions from ^spin^+i aj^fj ^spm.-i added. The contri- 
bution from jS'^p'"''^ vanishes anyway. Consequently, the 
result for the full integral is 



' d^P <^spin ^ SjS + 1) A^c 

N{0)T ^ N{0)S{S + 1) 

2Tspin ffcMB^i 6Tspin 



[go/J-BBe - giUBBi). 

(A8) 



Note that we have expressed this result in terms of /ii 
instead of the occupation fractions fm- This can be done 
in all our results so that a closed set of equations for the 
two magnetizations /ie and /ij is obtained. 

The integrals required for the equation of motion of /ic 
are quite similar. In the equation for the magnetization 
current we need integrals such as 



1 rd^pd^p' ^ 

~ N{0)t J (27r)6 ^ 
1 f d^p 



T J (27r)' 



(TVp Upa 



1 j/x 



T 5cAiB 



(A9) 



where the term with rip/g. vanishes since it is odd in p. 
Similar evaluations are required for and 5**^'". 



APPENDIX B: HYDRODYNAMIC EQUATIONS, 
VALENCE BAND 



Even though we restrict ourselves to the heavy-hole 
band, the angular integrals are much more complicated 
than in the conduction-band case since the explicit ex- 
pression H42(l for the hole magnetization /ih, the transi- 
tion probabilities, and the Zeeman energies now all de- 
pend on the direction in momentum space. As noted 
above, the analytical expressions for the transition prob- 
abilities are rather complicated. We use Mathematica to 
analytically perform the angular integrals of the form 



^cos"0|pO1A|/)p,|' 



(Bl) 



with n = 0,1,2 and A 
expressions like 

-'0\p{j\j^\j%, 



1, i"", j^, j , resulting in 



47r 



3 

40 



(7 - 6 cos 61' + cos 26*') cos^ — for / = =f3/2. 



{7 + 6 cos e' + cos 29') sin^ 



for / = ±3/2. 

(B2) 

Here, 9 and 9' are polar angles of p and p', respectively. 
As an example, we here evaluate the integral 

d^p 



f27r)3 



t^spm 



(B3) 



which corresponds to the one considered in App. ^ The 
collision integral is again divided into ^spm.o _j_^spin,+i _|_ 
^spin,-i rpj^g contribution from 5''P'"'° vanishes in anal- 
ogy with Eq. (jA3|) . The other terms are expanded in 
A/m up to linear order, 



■spin.itl 
pjm 



5](0) 



(B4) 



In the delta function is expanded in Bh, Bi and the 

noflip ^mp' where in the first 



term is then divided into y!-^'^ 



(second) term j' = j {j' — —j). The first term is eval- 
uated similarly to the conduction-band case, taking the 
more complicated angular integrals (jBip into account. 



.(0) 
"'nofiip 



S{S + l) 



2 ffhMB 



ffhMB^h 



-iV(0)5iMBi?i 

.(0) 



(B5) 



Also, writing out Sjjjp and renaming j <^ —j in the first 

E'^^ can also be 



term one can see that EL 



(0) 



.(0) 



flip noflip 

evaluated similarly to the conduct ion- band case, 



N{0)T 5 ^ 

'-'^'spm ^ 

^ m 



[-m±3m^TS[S+l)\ A/„, 



(B6) 



15 



which simphfies under summation over the three contri- 
butions, 



(2^)3 

N{Q)S{S + 1) 
^ 108t«, 



^(^+1) Mh 



18r, 



spin 5hA*B 

5A^(0)T Hi 



' spin 

m{Q)s{s + i) 



36r, 



spin ffiA^B'T'i 



108t, 



(B7) 



spin 



In the integrals pertaining to the hole magnetization and 
magnetization current we obtain some terms in which the 
occupation fractions fm cannot be reduced to /im- These 
terms cancel in the final equations of motion so that again 
a closed set of equations for /ih and /Xj is obtained. 



APPENDIX C: ABSENCE OF BERRY-PHASE 
CONTRIBUTIONS 



it commutes with the kinetic energy in the absence of 
spin-orbit coupling. Thus all terms in Eqs. HC1|I . HC2|I . 
and ljC3p vanish. 

For the valence band in the spherical approximation, 
Eq. H31|l . the spatial part of the Bloch wave function is 
also constant but the spin part is not. It is given by 
Eq. The p gradient is then 



Vp|j>p 



4> 



p sm f 



(C4) 



Furthermore, the Zeeman term does not commute with 
the kinetic energy so that we expect contributions from 
the perturbation. We use a perturbation expansion in the 
effective field to obtain the terms appearing in Eqs. (|C1(I 
and ljC2p . The hole Hamiltonian in the spherical approx- 
imation reads 



In the present appendix we show that Berry-phase cor- 
rections do no contribute to the hydrodynamic equations 
to linear order. In the framework of semiclassical the- 
ory they have been discussed in detail by Sundaram and 
Niu.-^ If one considers a wave packet made up of electrons 
of a single band, with narrow spread in real and momen- 
tum space, and with center-of-mass position r and mean 
momentum p, then the semiclassical equations of motion 
for these quantities are, in the absence of scattering)^^ 

r = VpEp, - zp„ ((VpulV^u) - (V>|Vpu)) 
((VpM|V» - (V>|VpU» 
-i{{Vpu\dtu)~{dtu\Vpu)), (CI) 
P = -V,Ep,+ipo,{{Vru\V^u)-{V^u\V,u)) 

+ lf„((Vr^/|V?u)-(V>|VrM)) 

+ i{{VrU\dtu)-{dtu\VrU)). (C2) 

Summation over a = 1,2,3 is implied. \u) = \upa) is the 
periodic part of the Bloch wave function and Ep^r is the 
wave-packet energy, which also contains a Berry-phase 
correction^ 

Epa - Ep, - Im (VrUp^l • {Ep„ - VpUp,), (C3) 

where He is the local Hamiltonian for the wave-packet 
center and momentum and Ep^ is the corresponding 
eigenenergy. This expression applies to conduction-band 
electrons; in the hole case a should be replaced by j. Note 
that the spatial gradient Vr acts on the center-of-mass 
vector, on which the states \u) depend parametrically. 

For the conduction band we can immediately see that 
Berry-phase effects are absent: In field- free equilibrium 
all spatial and temporal derivatives vanish. The p gra- 
dients also vanish since for the Hamiltonian H'^^^ = 
p^/(2mcb) the periodic part \u^^^) of the Bloch wave func- 
tion is constant and the spin part |±l/2) is also indepen- 
dent of p. This is not changed by the Zeeman term since 



i [(71 + 572/2) p2- 272(p-j)2] +5h^BS-zBh, 



Hi 



The unperturbed eigenenergies are 



(C5) 



2m 



P3 



(71 - 272) 
(71 + 272) 



(C6) 



for j = ±3/2, 
for j = ±1/2 

where only the spin part 
has a nontrivial p dependence. Assuming an effective 
field in the z direction, the first-order perturbation is 



and the eigenstates are ju'"-' 



P3 ■ 



(1) 1 

Cpj = 5hMB 3 



Restricted to heavy holes, e'^j = gh/^B (i/3) cosOBh. 
Degenerate perturbation theory yields the perturbations 
to the states, 



p,±3/2/ 



9h^J■B g -Bh 



\"p4/2l-' l"p,±3/2/ I (0) V 
(0) _ (0) l%,l/2/ 

''p,±3/2 p,l/2 



,(0) 



-1/2 



\u 



(0) N 
p,±3/2/ 



JO) 



JO) 



^p,±3/2 ^p,-l/2 



-1/2/ 



(C8) 



Introducing the difference between heavy- and light-hole 
energies, ^p = — 272p^/m, we obtain 



(1) > _ ffhA^B sin6'i3h , (0) , 
^P.±3/2) " 2V3.gp 



(C9) 



Simplifying the notation by writing only the spin part of 
the wave function, this gives 



p,±3/2/ 



ghHB sm9Bh 



2V3 



9p 



^2, 



(CIO) 
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The Berry-phase correction for the energy of heavy holes, 
given in Eq. HC3() . is, to first order. 



,(o)\ 



= Im 



,(0) 



P./3l(4'^^-^0)|Vp<,.,. 

(Cll) 



,(o)\ 



2V3 



Using that l^p^^yg) is an eigenstate of Hq we obtain 



9p 



ASpj = ^ • im 



4> %/3 



psm 



e 2 



smf 



P 

4p 



%/3 



(z X p) • VrSh- 



(C12) 



This correction evidently diverges for small p. The origin 
is the breakdown of perturbation theory as the energy 
difference gp between heavy and light holes goes to zero. 
This divergence is not crucial here since states deep inside 
the Fermi sea do not contribute to the response. 

The energy entering the semiclassical equations of mo- 
tion is, to first order, Epj = e^^ 
Eq. I|C2|I reads, to first order, 



.(1) 



Aepj . Thus 



= -gh^J■B 2 C0S6' Vr-Bh 



Vr [(Z X p) • VrBh] ■ 

(C13) 



Thus we find an additional force which is proportional 
to a second derivative of the field but independent of 
the spin direction, i.e., an orbital contribution. Then 
Eq. HC1|I becomes, dropping subscripts p, j, 

r - Vp^-zp„((Vpj.W|V>(°))-(V>W|Vpu("))) 
-iV^e(°) ((Vpu("'|V» - (V>|Vpu(°))) 
- i ((VpM(°) \dtu) - (9,^.1 VpM(°))) . (CM) 

The term multiplying pa can be evaluated explicitly and 
is found to vanish for the heavy holes. Thus 

• - P • D ^ 

r = .ghMB 7: sm - 

mhh 3 p 



ImVp(VrZ.|-(6(°)-i/(°))|Vp./°) 



u\VpU^"^)) 



- i ((Vpw(°) - (ft w| VpwW)) , (C15) 
cf. Eq. (|35|l . To first order, the Boltzmann equation reads 

dtnp, + -^-Vrnp,+p-Vpn^^^^S^f+J2s;f:}, (C16) 

since Vr/ipj and p are both linear in the perturbation. 
Thus the correction terms in Eq. ljC15|) drop out here and 



the only new term on the left-hand side comes from the 
orbital force in Eq. HC13(I . The equation of motion of /ih 
is obtained by multiplying the Boltzmann equation with 
~5h/^B(i/3) cos and summing over p, j. The orbital- 
force term drops out since it contains j — 0. 

The right-hand side of Eq. (|C16|I also has to be mul- 
tiplied with — (7hA'B(j/3) cos 6* and summed over p, j. 
The Berry-phase correction Aep^ to the energy appears 
in the delta functions implementing energy conserva- 
tion. If we evaluate the resulting integrals by expand- 
ing this delta function as in App. El all terms multi- 
plied with Acpj should be evaluated to order zero. Then 
the only j, j' dependence comes from the explicit fac- 
tor j/3 and from the transition probabilities. However, 
explicit evaluation in the 4x4 spin space shows that 
E,y073)lp(j|/)pf = 0, E,y(j73)lp(jb1/)pf = 0, 

E,y073)(lp(j|i+l/)p'P + lp(j|j-|j'>P'P) = so that 
all these terms vanish. Thus there is no contribution to 
the equation of motion for the hole magnetization. 

The equation of motion for the magnetization current 
contains an additional factor of r in the integrand, 

which should be calculated to linear order, see Eq. (jClSp . 

For the left-hand side we obtain 



.gh^B J (27r) 



d'^p 



ffh/^B I cos 9 VrSh • VpTl^"^ 



ffhMB 

4p 



Vr [p • (Vr X BiO] • Vpn<^^ 



(C17) 



Since all terms multiplied by r are already of first order 
we replace r by p/mjih. Then the first two terms in the 
parentheses are identical to the ones calculated above 
and the third vanishes due to J2j 3 — 0- Thus the left- 
hand side of the equation of motion is not changed by 
Berry-phase contributions. On the right-hand side we 
have to multiply the collision integrals by (j/3)cos6'r 
with r ^ p/mhh + Av from Eq. (jC15|l . Av contains the 
term from the p dependence of the Zeeman energy as well 
as the Berry-phase corrections. The contribution from 
p/mhh is what we have calculated in Sec. Ill Bl except for 
the additional Berry-phase correction Aepj in the delta 
functions. This correction is irrelevant, however, by the 
argument of the previous paragraph. 

In the second contribution, Av is of first order so that 
the collision integrals should be evaluated to order zero. 
But these are of course zero since there is no net scat- 
tering in equilibrium. Consequently, the linear contri- 
butions to the velocity r, in particular the Berry-phase 
corrections, drop out of the equation of motion for the 
magnetization current . In conclusion, we have shown 
that the hydrodynamic equations for the valence-band 
case are unaffected by Berry phases to linear order. The 
results of Sec. HI Bl) are thus correct. 
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